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Abstract 

I— I We derive a model for the propagation of short pulses in nonlinear me- 

^0 dia. The model is a higher order regularization of the short pulse equation 

Ph (SPE). The regularization term arises as the next term in the expansion 

pH of the susceptibility in derivation of the SPE. Without the regularization 

• T-H term there do not exist traveling pulses in the class of piecewise smooth 

(— I functions with one discontinuity. However, when the regularization term 

1 I is added we show, for a particular parameter regime, that the equation 

supports smooth traveling waves which have structure similar to solitary 
" ' waves of the modified KdV equation. The existence of such traveling 

pulses is proved via the Fenichel theory for singularly perturbed systems 
and a Melnikov type transversality calculation. Corresponding statements 
for the Ostrovsky equations are also included. 

cn 

O 1 Introduction 

00 

The propagation of pulses in a dielectric medium is normally modeled by the 
^ nonlinear Schrodinger equation (NLSE). There are two main assumptions in 

the derivation of the NLSE. The first assumption is that the response of the 
^ material to electromagnetic excitation attains a quasi-steady state, the second 

5^ assumption being that the pulse width is large in comparison to the scale of 

oscillation of the carrier frequency |21j . The experimental generation of shorter 
and shorter pulses means that the second assumption in the derivation of the 
NLSE may not be satisfied. Indeed it is possible to generate pulses whose width 
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is just a few cycles of the carrier frequency [TB]. Therefore, new models are 
needed for short pulse propagation. One approach to the propagation of ultra- 
short optical pulses is that taken by T. Schafer and C.E. Wayne in [27]. There 
the authors use the full one dimensional Maxwell equations and formally derive a 
reduced model by assuming an experimentally determined optical susceptibility 
coupled with an ansatz that captures short pulses. The equation they derive is 
the short-pulse equation (SPE) 

dtxu + au + d^u^ = (1.1) 

where u : {t, x) G M+ x E ^ M is the real component of the electric field in the 
direction transverse to the direction of propagation and a e M is a constant. It 



was shown there that ( 1.1 1 does not support traveling waves that are smooth. 



They then propose that this might reflect the Maxwell equations not possessing 



traveling wave solutions in the short-pulse regime or that equation ( 1 . 1 1 is an 



accurate model only for short times, and the breakdown of solutions occurs after 



the time of validity of ( 1.1 1 as an approximation to the Maxwell equations. We 
consider here another approach, namely that such traveling pulses exist if higher 
order terms in the approximation to the susceptibility are accounted for. 

First, we strengthen the non-existence result of |27j and show there are not 
even piecewise smooth traveling pulses to ( |1.1| with a single peak. The possi- 
bility of piecewise smooth traveling waves is motivated by the development of 
discontinuities for the Cauchy problem for ( |1.1[ ) . Typically, one cannot construct 



global smooth solutions to hyperbolic equations such as (1.1) since discontinu- 
ities in dtu and dxU tend to develop in finite time. Since traveling waves are 
globally-in-time defined objects, the development of these discontinuities is a 
definite obstruction to the existence of smooth traveling waves. Therefore it 
is natural to look for traveling waves whose profiles u{z) are smooth on either 
side of discontinuity occurring at some point z = z^, and at an algebraic 
condition is satisfied so as to ensure the piecewise smooth function is a distri- 
butional solution to ( |1.1[ ) . In section [s] we show that the SPE does not support 
traveling pulses in the space of piecewise smooth functions so that scenario (iii) 
above is also not possible. As explained in Section [3] the main obstruction to 
the existence of traveling waves is the discontinuity in the right-hand side of the 



profile equations ( 4.39[). Indeed, the discontinuity restricts homoclinic solutions 



to lie in the set — -y/c/S < m < +y/c/3, whereas the jump condition requires 
c = ul where c is the speed of the traveling wave and u■^, is the absolute value 
of u immediately to the right (or equivalently left) of the discontinuity. Clearly 
these two conditions are not compatible. We resolve this issue by proposing a 
regularization mechanism which has the effect of removing the discontinuity in 
the profile equations as well as the need for a jump condition. The nonexistence 
of traveling waves (whether smooth or with discontinuities) leads us to reinves- 
tigate the derivation of the SPE. By approximating the linear response function 
X^'^'' of the dielectric media to higher order we arrive at the regularized short 
pulse equation (RPSE) 

d^dtu + au + dlu^ + I3d^u = . (1.2) 
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which is the SPE ( 1 . 1 1 with a higher order dispersive term of strength p. 



We remark here that the RSPE is similar to the Ostrovsky equation 
[dtu + d^u^ + (idlu) + aw = 



(1.3) 



which was derived by L.A. Ostrovsky [23] as a model for internal solitary waves 
in the ocean with rotation effects of strength a < 0. The main differences being 
the cubic nonlinear term and both the signs and size of the physical parameters 
If (3 is small (1.2 1 is a singularly perturbed equation which induces a fast-slow 



structure in the equations that can be exploited by relating certain rescaled 
versions of the equation to other well known equations. Indeed, by setting 
T = t/ \/W\ Eind ^ = xj Vl/?!, one may rewrite ( |1.2| as 



dfdrU - 



sgn(/3)a|u = . 



(1.4) 



Setting /3 = in the RSPE with the usual scaling l\L^ yields the SPE while 
setting /3 = in the fast scaling (1.4) yields the modified Korteweg deVries 
equation (mKdV) 

dtu + da:U^ + sgii{l3)dlu ^ . (1.5) 

Hence in the slow scaling the RSPE is a small perturbation of the SPE while 
in the fast scaling the RSPE is a small perturbation of the mKdV equation. 
Given this, we expect that the traveling waves u{x — ct) of the RSPE, if they 
indeed do exist, are close to solutions of the SPE for |x — ct| ~ 0{1/ ^/\f3\), 
and close to traveling waves of mKdV for |a; — ct| ^ 0{y/\(3\) with a transition 
region for |a; — ct| ^ 0{^/\f3\). This scenario is proven via geometric singular 
perturbation theory in section |4.2[ While in the derivation of the RSPE (3 need 
not be small, one can nevertheless introduce a small parameter e (4.54) whereby 



under a suitable rescaling the equation is singularly perturbed in e. 



2 Derivation of the Regularized Short Pulse Equa- 
tion 

Consider the Maxwell equations in three space dimensions, 

^tB = -V X E, V • D = p ^ ^ 

^ 2.6 
9tD = VxH-j, V-B = ^ ' 

where E,H : {x,t) e M'^ x M.^ M'^ are the electric and magnetic fields, 
D, B : (x, t) g K'^ X i-^ M*^ are the electric and magnetic flux densities, p G K 
is the electric charge and j : {x, t) G M.^ x i-^ M.^ is the current density. To 



derive (1.2), we make several assumptions about the physical setup. The first is 



that the medium is a dielectric. This implies that the there are no free charges 



^We note that the regularized short pulse equation was proposed earlier in [22) in the 
context of plasma physics. In that context the equation make sense for any combination of 
parameters signs. 
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or currents. The second is that the dielectric medium is cubic and independent 
of space. These assumption are written as 

(HI) p = and j = 

(H2) B = /xoH and D EqE + eoP(a;, y, z, E) 
(H3) P(x, y, z, E) = x^'' * E + x^^^ * jEpE 

where /io, Sq are free space constants and * denotes convolution in time t. Hence, 
Maxwell equations become 



(a) 
ib) 



eodtE = V X H - eo^t {x^'^ * E + x^'^ * IE^e} 
HodtH = -V X E 



(2.7) 



We turn (2.7) into a wave equation for E, by combining the time derivative of 
([2j|(a), with (|2j|(b). Recalling the identity V x V x E = V(V • E) - V^E, 



we get 

V^E - V(V • E) ^ HQEodf (e + x^^^ * E + x^^^ * |E|^E 
We look for solutions of the form 

where : M x M+ M and E-^ : x ^ 1^3 jg g^^g^^ 



(2.8) 
(2.9) 



E^ 



{0,E2,E3) 



Substituting (2.9 1 into ([2.8^ yields 



= Moeoa? (x^^) * {v^^) + x"^^ * (t-3|E^pE^) 



(2.10) 

Since we are interested in the evolution of E along x, we can start by looking for 
solutions for which E-'-(?/, z, t) is constant vector. In this case, by simply redefin- 
ing the x^"*^ by constants we find that the equation for the transverse evolution 
of the electromagnetic field satisfies the scalar one dimensional equation 



a: 



(2.11) 



Here we view x as the evolution variable, so to put this in a form in which t is 
the evolution variable we make the change of coordinates {x, t) i— > (i, x) and set 
(MoEo)"^ Co to get 



(2.12) 
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where now * denotes convolution with respect to x. It is convenient to consider 
the Fourier transform of (2.121 which we write as 

c^ e 



^(3)„3 



(2.13) 



Remark 2.1 Since we are interested primarily in modeling the effect of the x^^^ 
term, for simplicity we will assume that the x^^^ response is instantaneous and 



modeled by the Dirac measure. Thus in (2.131 and what follows we set x*^'^^ — 1. 
We investigate the effect of different x^^^ response functions in future work. 



2.1 Approximating the x^^^ Term 

In dielectric materials, the x*^^-* term is modeled in Fourier space as [3} [21] 



2^, 



(2.14) 



Typically for silica fibers, there are three resonances of importance which occur 
at practical wavelengths, namely at wavelengths of A = 0.068^m, A = 0.116/xto 
and A = 9.896/zm. If one restricts attention to wavelengths between 0.25 and 



3.5 iim, approximate values for the various constants in (2.14) can be obtained 



by fitting experimental data for light propagation in silica to obtain [19] 



Xlilica(^) — 



0.696A2 



0.4079A2 



0.8974A2 



A2 - (0.0684)2 A2- (0.116)2 A2 - (9.896)2- 
where A is the wavelength 



A = 



277 
J' 



(2.15) 



(2.16) 



The idea in [27], developed further in this section, is that (2.151 can be 



expanded in powers of A, and that the expansion is in fact a good approximation 
over a particular range of wavelengths. Indeed, for wavelengths between 1.6/Lim 



and 3.0/im one can approximate (2.151 by 



(1) a.jl 
silica 



(A) ~ fl — a 



(2.17) 



for appropriate values of a and /i (see [27l |7| for a discussion on the validity of 
this approximation). Here we write an expansion of x''^^ beyond that of [27j by 
adding an additional term, 

2 



- (1) 
^sinca 



(A)=xira"(A)-:>(^)" 



IX — a 



(2.18) 
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For any fixed range of wavelengths, ( |2.18 1 leads to a better approximation to 
(2.151 since it contains three free parameters a,/i, instead of just two as in 
(2.171. Note also that since A = 2tt/£,, the additional i> term corresponds to a 



higher order derivative in physical space. Finally, due to the convolution term in 
the nonlinear Maxwell equation (and hence the product of x^^^ and u in (2.13 1), 



an implicit assumption in the derivation is that u is small outside the range for 
which x*-^-* is well approximated. 



Recalling ( 2.16 ) and plugging the expansion for x, 



(1) a,jl,i>f 



ilii 



and transforming back to physical variables leads to 



(0 (2.18) into (2.131 



dfV — j'^dlv — av + d^v^ + i^d^v 



where 



(2.19) 



a/cj 



(2.20) 



2.2 Parameter Values 

Up to this point we have made no mention of the signs or actual values of the 



parameters in (2.181. In this section we compute the signs and actual values 



of the parameters a, /2, and v which depend on the wavelength regime over 
which the x'^'' response function (2.15) is approximated. We follow the work of 



Schafer and Wayne [27], in making a least-squares fit to (2.15) over the range 



1.6/im < A < 3.0/iTO, see Fig 2.1. The sign of the parameter v in (2.18) leads 



to a parameter regime for which our results hold. This approximation does not 
capture the behavior of the response function over a larger region, in particular 
one that includes its upwardly convex part when A < 1. A least-squares fit 
over the wavelength range 0.25/xm < A < 3.5/iTO gives a value of v with the 
opposite sign. Our analysis does not cover this case, but we should point out 
that, although the approximation is better over the larger frequency range, it is 
not as good as the one in [27J, which we use here, near the carrier frequency of 
interest. 



2.2.1 Approximation over the regime 1.6/im < A < 3.0/zm 



Here we compute the parameters a,jl,v for the ansatz (2.18) by curve fitting 



over the exact same wavelength range as in the original paper [27] . The idea is 
that the extra free parameter v leads to a better fit, and in addition provides a 
non-zero higher order regularization term. In this case the parameters obtained 
by a least-squares fit are 



& = 4.287786709 x 10"\ jl = 1.112628599, v = 1.192774270 x 10"'' . (2.21) 



In light of (2.20) this leads to an RSPE equation (1.2 1 with a > 0,/3 > 0. 
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1.12 




2.5 
lambda 



Figure 1: Approximation of the response function for silica using (2.181 with 
a = 4.287786709 x lO^i, jl = 1.112628599, i> = 1.192774270 x 10^ The 



dashed hue is (2.18), the solid line is (2.15) 



2.3 The Short Pulse Ansatz 



To get (1.1 1 and (1.2) from (2.19) we introduce the short pulse expansion 
v{x,t) = Svo{x\t^)+6'^vi{x^,t^) + 

where 



X ~ jt 



-t. 



Let 



then substitution of (2.22) into (2.19) yields 

dtxVo + avo + dlvl + Pd'^VQ = o{S) 



(2.22) 
(2.23) 



(2.24) 



(2.25) 

where we have dropped the S superscripts on x,t. Then to solve (2.25) to order 
dt.vo + avo + dlvl + /39>o = . (2.26) 



S we need to solve 
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3 General Properties of the RPSE 

3.1 Conserved Quantities and Hamiltonian Formulation 

While both the SPE the mKdV equation are integrable [U [3S] so associated 
with them are an infinite number of conserved quantities, the arguments in 
[TUl [S] show that the a term in the RSPE destroys the integrable structure. 



Nevertheless there are conserved quantities associated with (1.2 1. These are the 
zero-mass law 



the momentum 



m{u):= / u{x,t)dx^O, (3.27) 
Jr 

M{u) := [ V? dx (3.28) 



and the energy 

H{u) ( - \ lul" + /3|9^up dx = const. (3.29) 

2 Jr 2 

With this definition for H , one can regard H as the Hamiltonian for the RSPE 
since one can write it as the gradient flow 

dtu = JVH{u) (3.30) 

where is the skew symmetric operator 

J = -a, . (3.31) 

3.2 Dispersive Regimes 

Consider the linear RSPE equation obtained by omitting the nonlinear term, 

dt^u + au + Pd^u = . (3.32) 



Plane waves of the form u — Ae'^^^ are solution to (3.32) provided that the 
pair {^,Lo) satisfy the dispersion relation 

c{0 = ^+Pe (3.33) 



where c(^) = w/^ is the phase speed. The dispersion relation (3.331 shows us 
there are two difi'erent dispersion mechanisms. For the sake of nomenclature, 
we call the a/^"^ term lower-order dispersion and the /3^^ term higher-order 
dispersion. There are essentially four dispersive regimes one can investigate. 
When both lower-order and higher-order dispersion mechanisms are absent (a = 
(3=0) the RSPE reduces to a Burgers equation with nonconvex flux which 
is known to develop shocks (c.f. [9]). When only the lower-order dispersion 
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mechanism is absent, the RSPE reduces to the modified KdV equation which 
has soHtary wave solutions (corresponding to homochnic orbits) 




z = x-ct, c,l3>0 (3.34) 



and travehng front solutions (corresponding to hetcroclinic orbits) 



= ±Vctanh l^y'— , z = x - ct, c>0,/3<0 (3.35) 

When only the higher order dispersion mechanism is absent, the RSPE equation 
reduces to the SPE equation. Finally, we investigate the regime where both 
dispersion mechanisms play a role, which is the RSPE equation. 



4 Solitary Waves of the SPE and RSPE 



4.1 Nonexistence of Traveling Waves for the SPE 

Here we prove that there do not exist piecewise smooth traveling pulse solu- 
tions to the SPE. By "pulse" here we mean that there is only one point of 



discontinuity. Rewrite ( 1 . 1 1 as the system 

dxq = u 
aq 



(a) 
(b) 



dfU + ag + dxU — 



(4.36) 



which can be thought of as an integrated version of (1.1 1 in the sense that one 
can obtain (1.1 1 by taking the spatial derivative of (4.36)(b). In this sense the 
quantity q would represent the electric potential of u. Written this way the 
hyperbolic nature of the SPE is apparent in that one can view it as a correction 
to Burgers equation with nonconvex flux. Given this, we do not expect the 



existence of smooth traveling waves for ( 1 . 1 1 since in general the initial value 



problem for equations of this type tend to form discontinuities in the derivatives 
dtu and dxU, which means we must look for distributional solutions. This has 
already been proven in f27'. Therefore, we look for traveling waves with profiles 
(g, u) that are smooth on either size of a point of discontinuity which are dis- 
tributional solutions to the original PDE The requirement that the trav- 
eling wave be a distributional solution yields an algebraic condition that must 
be satisfied at the discontinuity, which are the usual Rankine-Hugoniot jump 
conditions. Traveling waves with these types of discontinuities are phenomena 
in both hyperbolic evolution equations (c.f. [HEO]) and parabolic equations (c.f. 

mnaEQ]). 

It is not hard to see that piecewise smooth functions (g, u) with discon- 
tinuities across a smooth curve S :— {{t,x) : x — <j{t)} with unit normal 
u = (v^,!/^) = (1 -I- iT^)~^/^(l, — (t) are weak solutions of (4.361 if and only if 
( 4.36[ ) is satisfied in the classical sense on either side of S and the traces of (g, u) 
satisfy the jump condition 
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(a) M] = 

(6) a[[u]] = [[u 

where [[ • ]] denotes the jump across the curve S. 



311 (^-^^^ 



Note that (4.371 imphes q is continuous across S. In the case of symmetric 
Jumps, we can make an expHcit calculation. Let u^, ^ and suppose u(5 + 0) = 
It* = —u{S — 0). Then 

cr = > 0. (4.38) 

Suppose the pair {q,u){z) is smooth on either size of z = with nonequal 
left- and right-hand limits at z = 0. Then (g, u)(z) is a distributional solution 



to (4.36) if and only if the following hold 



1. (g,u) satisfy the profile equations 



q ^ u 

aq 
c — 3u^ 

pointwise in {z < 0} and {z > 0} 
2. The Rankine-Hugoniot conditions 

M] = 
c[M] = 

hold on z = 0. 



(4.39) 



(4.40) 



Now write the profile equations (4.391 and transmission condition (4.40) as 



w = F{'w) ± z > 

c[[w]] ^[[G{w)]] z = 



where 



F{w) := 



u 

aq 



(4.41) 



(4.42) 



The origin e MMs the unique fixed point of ( |4.41[ ). The linearization of F 
about the origin yields the linear system 



with 



A 



= Aw 

1 

a/c 



(4.43) 
(4.44) 
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so that it is a hyperbolic fixed point if 

a/c > 0. 



(4.45) 



Since G M is a hyperbolic fixed point of (4.411, there exist stable and 
unstable manifolds W{0), >V''(0) which locally can be described as graphs over 
the stable and unstable eigenspaces E"^ of the linearized system. These man- 
ifolds can be characterized the Hamiltonian structure of the profile equations. 
Define H{q, u) by 



Let 



(a) 
ib) 



H{q, u) — aq^ — i 

i;_ = {{q,u) e ir2 

E+ = {(<Z,ii) gM^ , 



< ^Jcji.qu < 0} 
|u| < \/^,qu > 0} 



and define the sets , i?+ by 

(a) 
(b) 



:={i/(g,u) = 0}nE_ 
H+ ■.= {H{q,u) = 0}n^+ 



(4.46) 



(4.47) 



(4.48) 



It is not hard to see that = >V"(0) and H+ = yy"(0). 

Next suppose that w(z) is a piecewise smooth traveling wave. Then w(z) is 
defined through the stable and unstable manifolds (4.48) by 

w{z) = w^{~z)x{z<o} + w+{z)x{z>o} (4.49) 



where w+ and W- are solutions to the initial value problem (4.391 with initial 
data 



(4.50) 



and X is the characteristic function. We will show that w{z) defined by ( 4.49[ ) 
cannot simultaneously satisfy the jump conditions and at the same time lie 
within the stable and unstable manifolds for every z S M. By construction, w{z) 
is smooth away from z = and satisfies (4.41) on M/{0}, thus w{z) satisfies 
the profile equation (4.41) pointwise on IR/{0}. Notice that by the Rankine- 
Hugoniot jump conditions and the fact that {q^,u^) must lie on yV"*(0), w has 
symmetric jumps, therefore the speed c of the traveling wave w{z) is 



c = w: > 



(4.51) 



However, for to he in W{0) we need 



3>< 



(4.52) 



which contradicts (4.51). 
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4.2 Existence of Solitary Waves of the RSPE and Ostro- 
vsky Equation 

Loosely speaking, the mechanism that leads to the nonexistence of traveling 
waves for the SPE is the discontinuous right-hand side of the profile equations. 
As shown in Section |4.1| the main obstruction to the existence of piecewise 
smooth traveling waves to the SPE is that in order for the the jump condition 
to hold, the value of u at the jump must be larger than the allowable range 
|m| < \/c/3 dictated by the singularity in the right-hand side of the profile 
equations. If one could eliminate this discontinuity, then it might be possible 
to obtain traveling waves. In order to eliminate the discontinuous right hand 
side of the profile equations, we introduced a small higher order term into the 
equation. This higher order term arose in expanding the susceptibility further 
so as to eliminate the zero wavelength singularity inherent in the derivation of 



the SPE. Here we study the regularized short pulse equation (1.2 1 and show 
that higher order dispersion is needed to obtain smooth pulses. For the sake of 
generality we will prove the theorem for general pure power nonlinearities, 

dta;U + au + dlu^ + Pdlu = 0. (4.53) 

which includes both the Ostrovsky equation (k=2) and RSPE (k— 3). 

We now state the main theorem of this section, the proof of which will oc- 
cupy the next several subsections. 



Theorem 4.1 (Existence of Traveling Waves for the RSPE) Consider the 
RSPE with general pure power nonlinearity (4.531 and let 



(4.54) 



1. Suppose k is an even integer and e is given by (4.541. Then for e > small 
enough there exists a traveling wave solution 



provided that 



u = u{z), z = X — ct, u{z) ^ as z — > ±(X) 



sgn(c) = sgn(a) = sgn(/3). 



(4.55) 



(4.56) 



2. Suppose k is an odd integer and e is given by (4.54). Then for e > small 
enough there exists a traveling wave solution 



u — u{z), z = X — ct, u{z) as z — > ±oo 



(4.57) 



provided that 

sgn(c) = sgn(a) = sgn(/3) = +1. 



(4.58) 
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Remark 4.1 Note that when k is an odd integer, as in the RPSE, we are only 
able to cover the physically relevant case (3 > Q (see Section 2.2.1). The case 



P < is also physically relevant for RSPE (see Section 2.2, 0.25iJ,m < A < 
3.5fj,m) but it not considered herein. 

Remark 4.2 The small parameter e is a ratio of three parameters and can be 
thought of as being made small by fixing two of the parameters and varying one 
of them. Hence by this rescaling our analysis of traveling waves of the RSPE 
covers the three cases (i) a small, (ii) (3 small, and (Hi) c large. 

Remark 4.3 Similar results can be found in \17\ \18f . There the authors use 
variational arguments to prove the existence of a ground state traveling wave to 
the Ostrovsky equation for small speeds c. The GSPT approach used here forms 
the basis of work in progress JBf on the construction of multi-pulse and periodic 
traveling waves which are not ground states. 



4.3 The Scaled Profile Equations and Fenichel Theory 

Traveling wave solutions 

u{x,t) = u{z) z := X — at (4.59) 



to (4.53) must satisfy the singularly perturbed fourth order equation 



dz2 V ^~d^ + - cu I + au = 



For any a, /?, c 7^ consider the rescaling u(z) — u(z) where 




Then upon dropping tildes the profile equations (4.601 become 



dz"^ 



dz^ 



sgn(/3)u'^ — sgn(c/9)M I + sgn(a/3)u ~ 



(4.60) 



(4.61) 



(4.62) 



where e is defined by (4.541. In order to use the GSPT framework we must 



rewrite the fourth order singulary perturbed equation (4.621 in standard fast- 
slow form. In doing this, a subtle point is the correct choice of slow and fast 
variables, and a correct identification and placement of the small parameter in 
the system of equations. The rewriting of the equation as a system is not unique, 
and neither is the placement of the small parameter in the equations, and hence 
it follows that the identification of fast and slow variables is not unique either. 



With this in mind, let e be as in (4.541 and set 



u — u, V — eii, w — ev + sgn/Jw'^ — sgn(/3c)u, y — w (4.63) 



SOLITARY WAVES OF THE RSPE AND OSTROVSKY EQUATIONS 14 



and consider u, v the fast variables and w, y the slow variables. The motivation 



for this choice of fast and slow variables (4.631 will become clear in the subse- 



quent analysis. We may then write the profile equations for traveling waves of 



(4.531 as the equivalent problem 



eu 



ev = w + sgn(c/3)u — sgn(/3)u''' 
w = -y 
y = sgn(a/3)u 



(4.64) 

where • denotes differentiation with respect to the slow variable z. Define a fast 



variable ^ = z/e, then (4.641 becomes 



(4.65) 



u — V 

v' — w + sgn(c/3)M — sgn(/?)u'' 
w' ~ —ey 
y' — esgn(a/3)u 

where / denotes differentiation with respect to the fast variable — z/e. We 
call (4.641 the slow scaling and (4.651 the fast scaling. 



We prove the existence of homoclinic orbits to ( 4.64| ) (or equivalently (4.65)) 
via the Fenichel theory for singularly perturbed systems of ordinary differential 
equations, which forms the basis of geometric singular perturbation analysis. 
We now briefly review the Fenichel theory. For a very readable account of this 
theory and its applications we refer to the C.I.M.E. lectures by Jones [15] or 
Szmolyan [25]. 

Consider the system of autonomous ordinary differential equations written 
in standard fast-slow form 



ex 

y 



I{x-,y;£) 
9{x,y;e) 



(4.66) 



where x G M", y G K™, /, e C°° and < e < Eq < 1- Let z denote the indepen- 
dent variable in (4.661, which is referred to as the slow scale. By introducing a 
fast scale ^ = z/e in (4.661, one obtains the equivalent system 



X 

y' 



fix,y;e) 
eg{x, y; e) 



(4.67) 



where ' = The basic idea of GSPT is to analyze (4.661 by combining 

information derived from the reduced problem 



and the layer problem 





y 

x' 

y' 



f{^,y;0) 
9{x,y;0) 

f{x,y;0) 




(4.68) 



(4.69) 
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which are the formal e = hmiting problems of (4.661 and (4.671. The fun 



damental connections between the formal limiting problems and (4.661 were 
laid down by Fenichel in [llj and consist of three main theorems referred 
to as Fenichel's First, Second, and Third theorems. The basic idea of the 



First Theorem is that (4.68) may be seen as a dynamical system on the set 
Ai — {{x,y) G M" X M™ : f{x,y;0) = 0} which is an invariant manifold of 
fixed points for (4.69). Because of the trivial dynamics of y' — 0, A4 cannot 



be a hyperbolic set with respect to the full flow of ( 4.69[ ). It can however be 
hyperboUc with respect to just the flow of x' — f{x,y; 0), which leads to the 
concept of normal hyperbolicity. 



Definition 4.1 (Normal Hyperbolicity) A manifold A4 is said to normally 



hyperbolic with respect to a flow of the fast-slow system (4.66) if the linearization 
of (4.66) at any point p Cz A4 when evaluated at e = has exactly m zero 



eigenvalues and n eigenvalues /i with 5R/i ^ 0. 

It turns out that this assumption is enough to prove the following persistence 
theorem of N. Fenichel, which guarantees that M persists as a manifold for 
e > small enough. 



Theorem 4.2 (Fenichel's First Theorem jllj . |15| ) Suppose A4 is normally 
hyperbolic. Then for e > but sufficiently small, there exists a manifold Ai^ 
that lies within 0{e) of Ai and is diffeomorphic to M.. Moreover, it is lo- 
cally invariant under the flow of (4.66), (4.67), and , including in e for any 
k < +00. 

Attached to each point p E M there is a one dimensional stable manifold 
>V*(p) and a one dimensional unstable manifold 'W''{p). We collect these man- 
ifolds together to form stable and unstable manifolds for the full limiting slow 
manifold M, given by 



= UpeMW"(p) 



(4.70) 



To address the question as to whether these structures persist when e > we 
turn to Fenichel's Second Theorem. 



Theorem 4.3 (Fenichel's Second Theorem 115] ^ Suppose A4 is nor- 
mally hyperbolic. Then for e > sufficiently small, there exist a manifold 
>V*(A^e) (resp. 'W'^{A4s)) that lie within 0{e) of, and is diffeomorphic to, 
^(A^) (resp. W"{A4)). Moreover, it is locally invariant under the flow of of 
( 4.66[ ), (4.67) and , including in e for any k < +oo. 



Moreover, there exist stable and unstable invariant foliations with base Ai^ 
with the dynamics along each foliation being a small perturbation of a suitable 
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restriction of the dynamics of (4.691 



Theorem 4.4 (Fenichel's Third Theorem [11', '15]) Suppose Ai is normally 
hyperbolic. Then for e > sufficiently small there exist for each point Pe G A^e, 



manifolds that form an invariant family relative to (4.661, (4.671, which we 



name W^{p) and W"(p), which are within 0{e) of and difjeomorphic to the 
corresponding manifolds when e = 0. 



4.4 Construction and Analysis of the Solitary Waves 
4.4.1 The Slow Manifold and Reduced Problem 



When e = in the slow scahng (4.641 we obtain 






w 

y 



w + sgn(c/?)M — sgn(/3)u'^ 

~y 

sgn{aP)u. 



(4.71) 



Thus there is a manifold of fixed points for (4.711 which we call the critical 
manifold Aio given by 



Mo = {{u, v,w,y) eR^ X I xR : tp{w),v = 0} 
where 'ip{w) is implicitly defined by the algebraic equation 
w + sgn(c/3)i/' - sgn{/3)tl;^ = 

and 

I := {u : sg-n{P)ku''~^) < sgn(c/3)} 



(4.72) 



(4.73) 



(4.74) 



This means that the flow on the limiting slow manifold Mo is given by 



w 
V 



sgn{aP)'ip{w). 



(4.75) 



Proposition 4.1 (Regime of Normal Hyperbolicity of A^o) The critical 
manifold Mo defined through (4.721 is normally hyperbolic with respect to the 



flow of (4.651 



Proof. Setting e = in the linearization of (4.651 at any point po £ Mq yields 
the matrix 

/ 10 0^ 

sgn(c/3) - sgn(/3)/cu'=-i 1 


\ / 



(4.76) 
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which has characteristic polynomial 

fi^ (/i^ - sgn(c/3) + sgn(/3)fcu'=-i) = 0. 

Hence zero is an eigenvalue of multiplicity two, and there are two distinct eigen- 
values 

M± = ^sgn(/3)(sgn(c) - ku''-'^) 
which have nonzero real part provided that 



sgn{P)ku < sgn(c/3). 



(4.77) 



However, by (4.741 every point in Mo satisfies (4.771 



□ 



Remark 4.4 Since u is multivalued as a graph over w, the interval I (4.74) is 
introduced to choose the branch which contains the origin. 



Proposition 4.2 (Existence and Characterization of Me) Under the hy- 

(4.78) 



pothesis of Proposition (4.1) there exists a manifold M^ given by 

Me = Mo + 0{e^) 



and the flow on Me is given by 



w = -y 

y = sgn{aP)ip{w) + 0{e'^) 



(4.79) 



Proof. For the regime of normal hypcrbolicity of Mq, Theorem |4?2] guaran- 
tees the existence of a manifold Me which is an 0{e) perturbation of A^o- A 
straightforward calculation shows that Al^ is in fact an 0{e'^) perturbation of 
Mo- □ 

Proposition 4.3 (Analysis of the Flow on Me) The origin of the system 



(4.71) is hyperbolic when viewed as a dynamical system on Me if and only if 

sgn(a) = sgn(c). 



(4.80) 



Furthermore, when (4.80) holds the origin has a one- dimensional stable and one- 



dimensional unstable manifold with eigenvectors of the linearization of (4.791 
at the origin begin O(e^) perturbations of 







/io = (double eigenvalue) 
( * \ 



■n- = 



/ 

v+ = 



\ sgn(Q:/3) / 



(4.81) 
(4.82) 



\ sgn(a/3) / 
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Proof. First note that the origin is an element of A^o- Since the origin is a 
fixed point for the full problem (4.641 it remains a fixed point for the flow on 
A^e also. To establish hypcrbolicity of the origin as well as the eigenvectors 
of the linearized flow, we need only establish them for the limiting problem 
( 4.71[ ) or (4.751 which are fj, = ±-y/— ■0'(O). By differentiating (4.731 we get 
ijj'iO) — ~sgn(a/3)/sgn(c/3), so that /x = ±\/sgn(Q;)sgn(c) and the result follows. 

□ 



4.4.2 Analysis of the Layer Problem 



We now consider the fast scaling (4.65 1. Setting e = we have the layer problem 

(4.83) 



u — V 

v' — w + sgn(c/3)u — sgn{(3)u'' 

w' ^ 

■■' = 



y 

Essential in the construction of the traveling wave is that the layer equation has 
an orbit that connects points on the critical manifold A^o- By the definition of 
Alo (4.721, this means the fixed points satisfy 

w ^ sgn(/3)M^ - sgn(c/3)-u*. (4.84) 

However, the fixed points of the layer problem must remain fixed points of the 
full system, and since we are looking for solutions homoclinic to the origin, we 
require w ~ which implies 



sgn(c). 



(4.85) 



Note that if k is odd the speed c must be positive. Thus (4.831 becomes 



u 



V 



v' — sgn(c/3)u — sgn(/3)M* 
w' = 
V' - 



(4.86) 



If /fc is odd the fixed points of (H^are (-1,0,0,0), (0,0,0,0), (1,0,0,0) while 
if k is even the fixed points of (4.86 1 are (0, 0, 0, 0) and (sgn(c), 0, 0, 0). 



Proposition 4.4 (Analysis of the Layer Problem) The layer problem (4.86) 
has an orbit homoclinic to the origin if and only if 



sgn(/?) = sgn(c) 



The orbit is given by 

MC) = ('Z(C),-'Z(C)tanh 



k-1 



C ,0,0) 



(4.87) 



(4.88) 
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wht 



q{0 = sgn(c) 



(fc + 1) 



sech'=-i 



fc-1 



-c 



(4.89) 



and sgn(c) = ±1 when k is even while sgn(c) — 1 when k is odd. 

Proof. Since the center directions play no role in the construction of the ho- 
moclinic orbit we can restrict the flow of (4.861 to the two-dimensional (u,w) 
phase space 



u 
v' 



sgn(c/?)M — sgn(/3)u'' 



(4.90) 



This equation is the profile of the generalized KdV (gKdV) solitary wave, which 
has the well known homoclinic solution (4.891 provided the sign of c and (3 are 
the same. □ 



Remark 4.5 Split the parameter space into P = P_UP+ where P± :— {(a, /3, c) € 

: sgn(a) = sgn(/3) = sgn(c) = ±1,£ ^ 1}. From (4.851 we see that for k odd, 
sgn(c) = 1 so that (4.871 implies sgn(/3) = 1 and (4.80) implies sgn(a) = 1. 
Thus for k odd we can construct solitary waves in the parameter space P_|_ while 
for k even we can construct solitary waves in V = P_ U P-|-. Without loss of 
generality, in what follows we assume (a, /3, c) G P 



4.5 Tangent Spaces and the Transversality Calculation 

Now that we have analyzed the e = limiting systems in both the slow and 
fast scaling, we prove the existence of solitary waves to the full e > problem. 
We prove this via a reversibility argument together with a transversality calcu- 
lation, which essentially entails showing that certain manifolds associated with 
the slow and fast orbits are transverse at e = and thus transverse for e > 
small enough. 



4.5.1 The Reversibility Argument 

Here we formulate and verify a condition that, together with the transversality 
calculation in Section 4. 5. 2| below, proves the existence of a homoclinic orbit of 
the profile equations for the RSPE equations. This condition essentially follows 
from reversability of the dynamical system. 



Proposition 4.5 (Condition for the Existence of a Homoclinic Orbit) 

Suppose r(z) : K — > is an orbit of (4.641 (or equivalently (4.651 ) which sat- 
isfies 



lim r(z) = 

z — ^ — oo 



(4.91) 
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and for some z ^ z* , 

nnr(z*) is nonempty. (4.92) 
Then T(z) is a homoclinic orbit. 

Proof. First notice without loss of generality the translational invariance allows 



us to set z* =0. Next, since the profile equations in both the slow scaling (4.641 



and the fast scaling (4.65) are invariant under the transformation 

zi-^—z {u,v,w,y) i-^ {u, —v,w, —y) 

there exists a reversibility operator TZ{z) such that 

n{u,v,w,y){z) = TZ{u,-v,w,-y){-z) 
with FixTZ the two dimensional plane 

n := {{u, v,w,y) eR^ ■.v=^0, y = 0}. 

One may take TZ as 



n 






V 



-1 



\ 





-1 / 



(4.93) 



(4.94) 



(4.95) 



(4.96) 



Recall that the origin of (4.64), (4.65) is a hyperbolic fixed point with two stable 



and two unstable eigenvalues. Clearly then there is a V{z) such that F satisfies 
(4.91). Suppose that in addition it satisfies (4.92) at the point z = 0. This 



means we have constructed {F(z) : — oo < z < 0}. By applying TZ to this 
portion of F we can construct {F(z) : < z < oo}. By reversibility we have 
that F(z) = 7^F(-z) so that 



lim F(z) = lim 7^F(-z) 

2— >00 Z — >-00 

= 7^ lim F(-z) 
= 



(4.97) 



which by definition means F(z) is a homoclinic orbit. □ 
By Proposition |4.5[ all that is left to prove in order to establish the existence 



orbit that satisfies (4.91) and (4.92) 



of a homoclinic orbit to (4.64), (4.65) is that for e > small enough there is an 



4.5.2 The Transversality Calculation 



Consider the profile equations in the fast scaling (4.65). When e = 0, (4.65) 



has a homoclinic orbit Q — (5,5', 0,0) which satisfies (4.91) and (4.92). Show- 
ing that this holds also for e > small enough amounts to showing that the 
conditions for the implicit function theorem hold. This in turn amounts to 
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showing that the evolution of the two dimensional unstable manifold under the 



flow of (4.65 1 when e — projected onto the orthogonal complement of 11, n-*-, 



IS nonzero. 



Consider the equations in the fast scaling (4.65 1 in differential form notation, 
du' = dv 

dv' = dw + {1 — sgn{c)ku''^^)du ■ ('4 981 

dw' = ~£dy ' 
dy' = edu. 



where we have used Remark |4.5| to simplify the equations. Setting e = in 
( |4.98| ) yields 

du' = dv 

dv' = dw + (1 — sgn{c)ku'^~^)du 
dw' = 
dy' = 



(4.99) 



Consider the two form 



dviO A dyiO e A2 



(4.100) 



Recall that the two dimensional critical manifold Mq is given by 

Mo ■■= {{u,v,w,y) eR^ X I xR:u = ifjiw), v = 0}. (4.101) 

By the Fenichel theory, the unstable (resp. stable) manifold of Ai^ is com- 
pletely foliated by smooth curves referred to as Fenichel fibers. Each Fenichel 
fiber intersects Mq at a unique point called the basepoint of the fiber. Thus, the 
foliation is a 2-parameter family of one-dimensional curves. The important fea- 
ture of these fibers is that points on a fiber correspond to initial conditions that 
asymptotically approach the orbit on as z ^ co (resp. z — > oo) that passes 
through the basepoint of that particular fiber. Let /"{po) denote an unstable 
fiber contained in yV^{Mo) which has basepoint po G Mo and /"(po) denote a 
stable fiber contained in W^iMo) which has basepoint po E Mo- The Fenichel 
theory enables us to identify lower-dimensional invariant manifolds within these 
stable and unstable manifolds. Let 70 C Mo be an orbit on the slow manifold 
Mo satisfying lim2_^_oo 7o = 0, then it has its own unstable manifold, denoted 
by yV"(7o), which is simply the union of all unstable fibers which have their 
basepoints lying on 70. We now proceeds with the calculation. 

Let rji be a vector tangent to the reduced problem which we take to be 

771 = (7/;'(u'),0,-l,sgn(a)) € TW" . (4.102) 
Since lim^^±oo w{C) = then by continuity 

lim V'(tf) = ^-'(0) = -1 (4.103) 
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and we can take 

7?i = (-l,0,l,sgn(c)). (4.104) 

By definition, rji is tangent to the reduced flow at z = —oo. Let <i>^ = 
((^f , (^3, (^1) denote the flow of (4.991. Then 70 :— <i>^ • 771 deflnes an or- 
bit in Aio. Notice that for the layer problem (4.991, for every C S M there is no 
flow for both w and y so that the third and fourth component rji is invariant 
under <i>^, that is ■ iji — —tpf ■ (1,0, l,sgn(Q;)). 



Next take a vector tangent to the e = homoclinic orbit h (4.. 



take as the vector field of the fiow in the fast scaling at e = 0, (4.83) 
ri2 — {v,w + u — sgn{c)u'',0, 0). 



which we 



(4.105) 



Clearly, (4.105) gives a vector tangent to the homoclinic orbit h for every z e M 
and furthermore, 772 G TW". We now compute the projection of the limiting 
e = fast flow onto 11 at z = when applied to the the tangent space of W". 
That is, we wish to compute 



(dwAdy)(77i,r/2)(0). 



(4.106) 



To do this notice that 



{dv A dy) (ryi,772)(C) - {dv' A dy){m, 1^2) {0 + {dv A dy')(m, '72)(C) 

= (1 - sgn(c)/c«'=-i)(du A dy)i7ji,7j2m + (dv A 0) (771, 772) (C) 

= {1 ~ sgn{c)ku''-'^){du /\ dy){iji,r]2){0 
= -(1 - sgn(c)fcu''"^)i; 

sgn(c)fcu'^^"^)u' . 



Thus 



{dvAdy)im,r]2M = 



(4.107) 







00 dC 




(dvAdy) {T]i,r]2){0dC 

- f (l-sgn(c)fc9'-i)9'dC 

- ((?(0) - sgn(c)(7(0)^) 



(4.108) 



sgn(c) 



1 



fc+ 1 



which is nonzero since by assumption a, c nonzero. Since W" and yV| intersect 
when £ = (4.108) shows that this intersection is transverse. Therefore, by the 
implicit function theorem, the manifolds still intersect for e > small enough. 
The intersection of W" and for nonzero e finishes the construction of the 
pulse. 

□ 



SOLITARY WAVES OF THE RSPE AND OSTROVSKY EQUATIONS 23 



4.6 The Melnikov Calculation, Homoclinic Breaking, and 
Asymptotic Decay of the Wave 

Here we want to prove some analytic and geometric properties of the wave for 
\z\ large. In particular we want to show that the homoclinic orbit that exists 
when e > enters (resp. exits) tangent to the weakly stable (resp. unstable) 
eigenvectors. The idea is to show that when e > is small but nonzero, the 
homoclinic orbit (4.88) of the e = layer problem breaks. That is, the one 
dimensional stable and unstable manifolds of the origin which intersect for the 
e = layer problem, fail to intersect when e > 0. Thus the homoclinic orbit 
of the RSPE equations cannot enter (resp. exit) tangent to the eigenvectors 
associated with the linearization of the layer problem at the origin. This will 
be proved via a Melnikov calculation. Using the fact that the profile equations 
for RSPE are reversible, this means that the solution must enter (resp. exit) 
the origin tangent to the eigenvectors associated with the linearization of the 
reduced problem at the origin. 



4.6.1 The Melnikov Integral Calculation 

Here we show that the homoclinic orbit that exists in the fast scaling when 
£ = breaks when e > 0. To set up for the calculation, write the traveling wave 



equations in the fast scaling (4.621 as 



u 
v' 
w' 

y' 



w + u — sgn(c)u'^ 

y 



(4.109) 



9 

~e u . 



Notice here that once again the issue of the correct placement of the parameter 
e is important. Set U = {u,v,w,y) and write (4.1091 as 

U'^F{U;e^) (4.110) 



with F defined by the right hand side of (4.1091 
in Proposition 



4.4 



When e 



we have shown 

that the equations posses a homoclinic orbit h = (q, g',0,0) 
(4.881. By reversibility of the profile equations for the layer problem (also set- 



ting e = in (4.109)) we see that the one dimensional stable and unstable 
manifolds for the layer problem intersect in the (u, v) plane along v = 0. One 
way to measure how much the stable and unstable manifolds of the the layer 
problem miss each other is to define the distance between these curves evaluated 
along V = which to first order is given by the Melnikov integral (4.117) which 
we now describe. 



Consider the variational equations obtained by linearizing about the e = 
homoclinic orbit h = (g, g',0,0) (4.88), given by 



V' = DuF{h;0)V 



(4.111) 
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which we write expHcitly as 



v = 



V 





1 — sgn(c)fcq'^~^ 





1 








\ 





1 













1 













/ 



V. 



(4.112) 



The adjoint variational equations 



(4.113) 



are given expHcitly by the system 



/ 





— (1 - sgn(c)fcg'^ ^) 








\ 




-1 





1 










-1 










V 








-1 








(4.114) 



where ^' = (-01, ^^2, V'3j V'4)- From (4.114) we see that -02 satisfies the equation 

(4.115) 



1 - sgn{c)kq'' V2 



so that -02 = q' solves (4.1151 since it yields the profile equations for the gKdV 
equation (4.901. Since V'l = —"021 V'3 = —"02 and -04 = —03, 



03 
V ^-4 / 



/ -q" \ 



Let 



where here 



/oo 
-00 



(4.116) 



(4.117) 



• ) denotes the vector inner product, and F{h; 0) is the vector 
field defined by setting e = in the right hand side of (4.109). Note that while 
DuF{U;0) is not of full rank and thus the origin of (4.109) is not hyperbolic, 
M defined by (4.117) nonetheless defines a Melnikov integral in the usual sense 
(c.f. |1]). Thus the homoclinic orbit breaks for e > if Af 7^ 0. We now evaluate 
the Melnikov integral (4.117). 



Clearly 



D,.F{U;e')\ = 








(4.118) 
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so we have 

/oo 
{D,. F{h{C);e%^^MO) 
'OO 

/oo 
-OO 




(4.119) 



Since q > for all C € R we have 

M^O (4.120) 

which proves the result. 



Remark 4.6 The geometric significance of the Melnikov calculation is that the 
homoclinic orbit of the full £ > problem cannot enter (resp. exit) the origin 
tangent to the strongly stable (resp. strongly unstable) eigenvectors. This sit- 
uation is evidence of an orbit-flip bifurcation in e, and allows us to construct 
m,ulti-bump traveling waves 

Proposition 4.6 (Asymptotic Decay of the Wave.) Assume e > fixed 

(4.121) 



and u{z) is a solitary wave of (4.53). Let 



Me = iMri = Mr = 



V2sgn(c)-2Vl-4e2. 



Then for any m > there exists an R > large enough so that 
for \z\ > R. 



(4.122) 



Proof. The Melnikov calculation coupled with the reversibility argument means 



that the homoclinic orbit of (4.64 1 (4.65 1 which we constructed for e > small 



enough cannot enter (resp. exit) tangent to the fast directions, and so must 
enter (resp. exit) tangent to the slow directions. The eigenvectors associated to 
the slow directions have magnitude /x^. □ 
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